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1. Introduction and motivation 

Among several features of QCD the chiral symmetry breaking is one of the most interesting 
ones. At low temperature this symmetry is spontaneously broken and the vacuum develops a quark 
anti-quark condensate, (qq) ^ 0. Massless Nambu-Goldstone (NG) bosons should appear in the 
spectrum of the massless theory. In real QCD, of course, the light quark masses explicitly break the 
chiral symmetry, giving a small mass to the NG bosons that have been identified as the 8 lightest 
mesons (pions, kaons, T]). Classically, the pattern of chiral symmetry breaking is the following (Nf 
being the number of light quarks) 1 : 

SU(N f )v x SU(N f ) A x U{l) v x U(l) A -> SU(N f ) v x U(l) v (1-1) 

that should actually give 9 NG bosons, but the ninth particle is absent in the spectrum. This apparent 
problem was solved noticing that the flavor- singlet axial U(\)a rotation is no more a symmetry at 
the quantum level (even in the massless limit); it is anomalous [1] due to the presence of instanton- 
like configurations. A direct effect of the anomaly is the large splitting in the mass of flavor-singlet 
and non-singlet pseudoscalar mesons (see Witten-Veneziano formula [2, 3]). 

While at zero temperature the physics is quite clear, at finite temperature still there is no defi- 
nite answer to the question if axial U ( 1 )a symmetry is restored or not. If it is restored, an interesting 
problem is to establish if this happens at the same critical temperature of chiral symmetry restora- 
tion. This would have relevant effects on the pattern of symmetry breaking and so on the critical 
exponents of the phase transition [4]. 

By semiclassical calculations of dilute instanton gas at very high temperature T S> T c , we ex- 
pect a strong suppression but not an exact restoration of the U(1)a symmetry. The most advanced 
lattice result in this context is the one by Vranas [5], where, using domain wall fermions, he con- 
cluded that just above chiral phase transition the U(1)a symmetry remains broken but by a very 
small amount. It is an open question what is the effect that may have to the order of the transition. 

Our aim is to study the fate of chiral and U ( 1 )a symmetry (and the mass of particles) at 
finite temperature around and above the phase transition using the fermionic action that retains 
the maximal amount of chirality on the lattice, i.e. the overlap formulation [6]. The JLQCD and 
TWQCD collaborations have performed large scale QCD simulations using the overlap action [7]. 
All simulations were done at zero temperature and we investigated the chiral behavior of spectra, 
low energy constants, chiral condensate and topological susceptibility [8, 9, 10, 11]. A non-zero 
topological susceptibility, indicating anomalous breaking of U ( 1 )a symmetry, was clearly observed 
in those simulations at T = and also the linear dependence with sea quark mass was obtained, as 
predicted by chiral perturbation theory. 

In order to simulate QCD using HMC with dynamical overlap fermions, fixing the topological 
sector was crucial because allowing for topology change would be extremely expensive [12]. To 
run a simulation at fixed topological charge Q it was introduced an irrelevant term in the action 
[12], that suppresses the occurrence of zero eigenvalues of the hermitian Wilson- Dirac operator 
that have to be crossed in order to change topological sector. Fixing topology, of course, creates 
a bias in physical results that must be corrected. A full theory describing the effects of working 

l U(l)v gives the conserved baryon number and SU {Nf)y is only softly broken by the small quark mass difference. 
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at fixed Q was developed [13]: the effects at zero temperature are understood, under control and 
0(1 /V). We will discuss with more details the subject in the following section. 

In order to obtain reliable results at finite temperature we need to check whether the same 
methods used at zero temperature to correct for fixed topology effects work even in this case. So, 
we started with some exploratory studies using quenched theory but fixing topology in order to 
compare with previous results in the literature. We measured the topological susceptibility and 
several meson correlators at finite temperature. We will report the results of these simulations and 
the preliminary results in full QCD to investigate U ( 1 )a restoration. 



2. Simulations and results 

Before discussing the results, let us briefly describe the methods used to measure correlators 
and the topological susceptibility at fixed topology. Detailed derivation of equations can be found 
in [13]. 

By using a saddle point expansion of the QCD partition function in a finite volume we can 
derive an expression for Zg, partition function at fixed topology (V is the 4- volume): 



1 



a gaussian distribution for topological charge that can be used to show that: 



Q 2 



i-^ + o( 1 



(2.1) 



<p(*)p(0)> = 1 y(v- Xt -^v) + ° (y ~ 3) - (Z2) 

which implies that the topological susceptibility can be extracted from a long range conelation of 
the topological charge density p(x). At first order in 1 /V we can ignore the contribution of the C4 
term, and check later the consistency of the assumption. 

Using the overlap operator we can define an object that has the same properties as the topo- 
logical charge density: 

p m (x) = mtr[Y 5 {D c + m)-l], (2.3) 

where (D c + m)~ l is the valence quark propagator constructed using the chirally symmetric over- 
lap operator. An alternative way is to consider the pseudoscalar isosinglet v\' correlator, whose 
disconnected part is equal to (p m (x)p m (0)) at large distances and couples only to the fast decaying 
T]' state, making it a better choice in order to estimate the long distance limit. 

We measured the connected and disconnected part of pseudoscalar correlators at long distance 
to extract % t [8]. Since we are working at finite temperature the correlators are measured and aver- 
aged over spatial directions. We reconstructed the correlators by using the first 50 eigenvectors of 
the overlap operator assuming low mode dominance. For example, the connected scalar correlator 
is given by 

C(x,y) N = Tr^ — rr— - (2.4) 

and similar expression of the pseudoscalar (just iXi — > —iXi). We checked that the saturation with 
50 eigenmodes is sufficiently accurate for the infrared behavior. 
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2.1 Pure gauge simulations 

By introducing the topology fixing term in pure gauge simulations we can check if, even at 
finite temperature, we can reconstruct topological susceptibility using the method described above 
by comparing with the literature. 

The setup is the following: Iwasaki action + topology fixing term at temperatures ranging from 
[0.8, 1.3]r c on two different volumes 16 3 x 6 and 24 3 x 6. The critical point was estimated to be at 
j8 c = 2.445 by inspecting the Polyakov loop. 

We first check whether the eigenvalue distribution behaves as expected. The typical distribu- 
tion is shown in Figure 1. We do not find any discrepancy with previous results (e.g. [14]). The 
presence of a peak for small eigenvalues in the high temperature side (T > T c ) was confirmed. In 
[14] these modes are associated with the presence of dilute gas of instantons-anti instantons. 

The most interesting result at this stage is the behavior of the topological susceptibility at 
finite temperature in comparison with the results of Gattringer et al. [15], shown in Figure 2. The 
asymptotic value for the disconnected correlator (see equation (2.2) ) was estimated using a joint 
fit of the connected and disconnected parts and assuming a double pole form for the last one in the 
quenched theory. The decay is dominated by pionic states for both of them. Then we use (2.2) to 
extract % t assuming that the C4 term is negligible. 

In [15] the topological susceptibility was measured using the index theorem by just counting 
the number of zero modes of an approximated version of the overlap operator, which is a clean 
definition without the ambiguities due to cooling. By also performing a simulation without the 
topology fixing term, we checked their results also with our exact overlap operator, finding no 
significant deviations within errors (cross symbol in Figure 2, temperature T /T c = 1.1). 

Below the transition temperature there is agreement between the two sets of data, showing 
that the method works very well at least until T c . Above the critical temperature our results are 
systematically lower than the reference ones. We are currently investigating the source of this 
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Figure 2: Comparison of topological susceptibility results with [15], lattice 24 3 x 6 



discrepancy which could be traced back in the assumptions leading to (2.1) or, the c\ term could 
be non negligible in this regime. We are trying to estimate this quantity but still we are getting too 
large error bars for the four-point spatial correlators. 

Another possibility to check the validity of (2.2) is to measure it on different topological sec- 
tors. Simulations are on the way at the time of writing. 

2.2 Full QCD simulations 

In full QCD with two flavors of dynamical overlap fermions we concentrated on the channels 
71,8,1]', a given by the correlators of the operators yjsTy, ^/ti/a, YJsY, ty^y, respectively. Here, T 
is the Pauli matrix in the flavor space. 

If chiral symmetry is restored we expect that the pairs (o,7t), (i]',8) become degenerate. 
Similarly, if the flavour-singlet axial symmetry is restored too at high temperatures we would have 
that all the channels become degenerate. The % and the i]' differ just by the disconnected part, 
which is essentially given by the near-zero modes. This observation implies that the U ( 1 )a breaking 
is driven by near-zero modes. By looking at the spatial correlators in those channels, we can check 
whether it happens. The problem of establishing if the axial symmetry is restored at the critical 
point is still open and has a relevance on the possible order of the phase transition [4]. 

Some details on the simulation follow. The algorithm is HMC, using Iwasaki action with 
topology fixing term and two flavors of sea quarks. The size of the lattice is 16 3 x 8. We choose 
N t = 8 to ensure that the configurations are smooth enough. Masses start from am = 0.05 down 
to am = 0.01 giving a pion mass of around 400 MeV at lowest temperature. The j3s were chosen 
to be in the temperature range T = [171,243] MeV, T/T c = [0.97,1.39] (assuming the critical 
temperature to be T c = 175 MeV). A comment is in order here: we couldn't estimate directly the 
transition temperature because it requires long history runs to measure susceptibilities. Anyway 
we checked that, by looking at eigenvalues of the Dirac operator (Figure 3), and the Polyakov loop, 
we simulated /3s just below and above the transition temperature. The topological sector is mainly 
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Figure 3: Eigenvalues density in full QCD simulations, lattice 16 3 x 8. Density around zero gives the chiral 
condensate by Banks-Casher relation. 



Q = but we have some simulations also at Q = 2. All configurations were generated and analyzed 
using the BlueGene/L installation at KEK. 

Some preliminary results are shown in Figure 4. At this stage we cannot say anything quan- 
titative since the volume seems still too small to extract masses. We can just check by inspection 
of the plots when the correlators become equal. We observe that the correlators start to become 
almost degenerate after the transition temperature and that as the sea quark mass is decreased this 
degeneracy is improved. We are collecting more data points to extrapolate this result toward the 
transition temperature. 



3. Conclusions 



We described our project on finite temperature QCD with overlap fermions. Simulating over- 
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Figure 4: Scalar and pseudoscalar spatial correlators at finite temperature. The estimate for the a, green 
area, correlator has still huge errors in comparison to the others. We see some degeneracy for the correlators 
at J3 = 2.30, T = 208 MeV, right panel. 
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lap fermions forces to fix the topology in our case. A method to extract physical results from fixed 
Q simulations was previosly developed for the zero temperature regime. To obtain reliable results 
we must check that we can apply the same method at finite temperature or find the necessary mod- 
ifications. We started the investigation analyzing the behavior of topological susceptibility in pure 
gauge theory at Q = 0. We found that our results differ at high temperature from the previous 
works. We are currently investigating the source of this discrepancy. A deep understanding of this 
problem is essential in the interpretation of full QCD results where we found restoration of axial 
symmetry at least from temperatures above 1.1T C . 
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